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Quantum simulations of spin systems could enable the solution of problems which otherwise require 
infeasible classical resources. Such a simulation may be implemented using a well- controlled system 
of effective spins, such as a two-dimensional lattice of locally interacting ions. We propose here a 
layered planar rf trap design that can be used to create arbitrary two-dimensional lattices of ions. 
The design also leads naturally to ease of microfabrication. As a first experimental demonstration, 
we confine ^^Sr^ ions in a mm-scale lattice trap and verify numerical models of the trap by measuring 
the motional frequencies. We also confine 440 nm diameter charged microspheres and observe ion- 
ion repulsion between ions in neighboring lattice sites. Our design, when scaled to smaller ion-ion 
distances, is appropriate for quantum simulation schemes, e.g. that of Porras and Cirac (PRL 92 
207901 (2004)). We note, however, that in practical realizations of the trap, an increase in the 
\^ secular frequency with decreasing ion spacing may make a coupling rate that is large relative to the 

decoherence rate in such a trap difficult to achieve. 
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4^ I. INTRODUCTION 

Oh 

Elucidating the low temperature properties of magnetic 
materials remains a challenge to computational physics. The 
^ combinatorial number of states in a many-spin system, as 
^^jwell as the quantum- mechanical nature of the materials, 
makes simulations difficult. These effects are strongly depen- 
^-H dent on the geometry of the spin system and the existence of 
^ defects in the spin lattice [U El [3] . Current tractable models 
of these systems rely on insights into the geometry that yield 
models with fewer degrees of freedom [J. 

fS| An alternative to classical numerical computation is quan- 
tum simulation, in which one well-controlled quantum sys- 
tem is used to simulate the properties of another. The type of 

00 simulation possible is currently limited by the controllability 
of the system and the number of degrees of freedom; a trade- 
^ off between the two has existed in every experimental imple- 

• ^ mentation of quantum simulation to date. Small molecule 

^ NMR [i 1 H [8] and experiments with a small number of 
^ trapped ions [9l [TOj [IT] have demonstrated the principle of 
universal quantum simulation for small systems. On the 
other hand, atoms in optical lattices have been used to sim- 
ulate large systems but with limited control [121 HSl HH US] • 
Theoretical proposals for quantum simulation of spin sys- 
tems have included the addition of controls to current optical 
lattice experiments [16", TT^, H^S] or the use of many trapped 
ions [19, 20, 21, 22, 23, 24 . Porras and Cirac show in Ref. 
[2T] a way to generate an effective antiferromagnetic 2D Ising 



*Both authors contributed equally to this work. 

'''Email: robclark@mit.edu 

^ Email : t ongyan@physics .harvard . edu 



Photoionization Laser 




Atomic Beam 



FIG. 1: Schematic of the lattice trap. An array of traps is 
produced by a single rf electrode with a regular array of holes, 
mounted above a grounded electrode. Ions will preferentially be 
loaded from a broad atomic beam at the intersection of the cooling 
and photoionization lasers. 



interaction in an array of ions. Such a simulation has recently 
been performed for the ferromagnetic case with two ions in a 
linear trap [11 . However, the observation of spin frustration 
requires a 2-D array of ions. One of the principal technical 
questions of quantum simulation research has been how to 
realize such a two-dimensional lattice of ions. A 2-D array 
of ions has been realized in a Penning trap by NIST [25j, 
but inconveniently the crystal rotates due to the crossed E 
and B fields. Arrays of electrons in individual Penning traps 
[26^ have also been proposed as one solution, although the 
schemes of quantum simulation mentioned above for trapped 
ions are not directly applicable. A recent proposal [27] pro- 
poses a scheme for doing quantum simulations in 2-D ar- 
rays of microtraps using localized electromagnetic fields and 
magnetic field gradients. This scheme does not use pushing 
forces from lasers, as in Ref. [211, eliminating errors due to 
sponataneous scattering and simplifying the optics needed 
considerably. 

In this paper, we present a layered planar rf electrode ge- 
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ometry that creates a 2-D ion lattice for quantum simulation 
of Hamiltonians in 2-D. The ion trap consists of a planar elec- 
trode with a regular array of holes, held at a radiofrequency 
(rf) potential, mounted above a grounded planar electrode. 
A single ion is trapped above each hole in the rf electrode 
(Fig. [T]). Ions will be preferentially loaded above the trap 
electrode at the intersection of the Doppler cooling and pho- 
toionization beams, allowing the user to write an arbitrary 2- 
D lattice structure. Also, in this Paul trap-based scheme, we 
avoid the large Zeeman shifts associated with Penning traps, 
as well as the rotation of the ion crystal. Understanding the 
basic properties of this trap is important not only for Porras 
and Cirac's proposal, but also for related methods such as 
Chiaverini and Lybarger's magnetic field gradient approach, 
and any ion trap quantum simulation protocol that relies on 
a fixed 2-D array of trapped ions. Our research here focuses 
on three trap properties needed for quantum simulation: the 
ability to stably confine ions, predictable trapping potentials, 
and measurable interactions of ions located in adjacent wells. 
Therefore, in this work, we ask the following questions: 1) 
How well can our design be used to trap an array of ions? 2) 
How well do numerical models of the trap match its observed 
properties? and 3) What ion-ion interactions between ions 
in neighboring wells can be observed? 

The paper is organized as follows. We first present a theo- 
retical model of the lattice trap (Section H), and then report 
on two demonstrations of the trap. In the first, we confine 
^^Sr+ ions and test numerical models of the trap by mea- 
suring the motional frequencies of the ions (Section III). We 
then use 440 nm-sized charged microspheres to measure ion- 
ion repulsion (Section IV). In Section V, we estimate the 
ion-ion spacing required in our trap to realize a quantum 
simulation. In Section VI, we summarize our results and 
discuss future work. 



II. LATTICE TRAP DESIGN 



imated by a slow secular motion superposed with a rapid 
oscillation, the micromotion, due to the oscillation in the 
potential Vrf. For Uq = 0, time-averaging the ion motion in 
the secular approximation (g <C 1) gives the following quasi- 
static pseudopotential which governs the secular motion: 
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Here ^{x) is the electrostatic potential when the drive volt- 
age V is applied to the ring electrode. 

The lattice trap can be thought of as a planar array of ring 
traps. Ions are confined in a 2-D lattice of potential wells. As 
discussed in Sec. I, this trap comprises two layers: an rf plate 
(extended ring electrode), and a grounded plane beneath it. 
At the center of each trap, the electric field associated with 
the electrostatic potential ^{x) is 0. Assuming approximate 
rotational symmetry in the plane of the trap, ^{x) has a 
multipole expansion: 



^{x) = V- 
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where r is the radial distance from the central axis of the 
lattice site and z is the distance along the central axis. The 
above expression is valid for infinite lattices, but for lattice 
traps containing many ions the potential will be correct near 
the center lattice site. The z = plane is defined such that 
it coincides with the point of null electric field. Eq. |2] defines 
two constants which depend on the trap geometry: ri, with 
dimension of length, and a, which is dimensionless. 
The pseudopotential is given by 
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Our lattice trap is an extension of the three-dimensional 
ring Paul trap [28 . Following this reference, we first review 
the theory of the ring trap. The ring electrode geometry is 
shown in Fig. [2] An alternating voltage of the form Vrf = 
V cos fit is applied to the ring electrode and the endcaps are 
grounded. The equations of motion for an ion in the ring trap 
are a set of Mathieu equations which have regions of stability 
depending on the dimensionless parameters a — ^2^1 and 

q = ^^2^2 ' Q and m are charge and mass of the ion, Uq 
is any DC voltage applied to both endcaps of the ring trap, 
and the distance ro is the distance from the trap center to 
the rf electrode (as shown in Fig. [2|. When Uq = and the 
system is in vacuum, the condition for stability is q < .908. 
We assume the trajectory of a trapped ion is well approx- 
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FIG. 2: Cross section of a ring trap. Ions are confined at the 
center of the trap. The ring is held at an alternating RF potential 
relative to the endcaps. Here ro = V^zq and the endcaps and ring 
electrode are hyperbolically shaped. 
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From the pseudopotential we define secular frequencies which 
characterize the curvature of the pseudopotential in the har- 
monic region: 
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where uOj. is the secular frequency in the plane of the trap and 
uOz is the secular frequency perpendicular to the plane of the 
trap. Note that ^ q so that uOz/^ gives a direct mea- 

sure of the stability of the confined ions. The micromotion 
approximation is best when uOz 

An additional grounded plate may be added above the ions 
to shield them from stray charges, and a static potential U 
may be applied to it. This change can be modeled by adding 
an extra U{z — zo)/zi in the pseudopotential, where Zi is a 
geometric factor with dimensions of length that depends on 
the height of this plane above the rf electrode and is com- 
puted, in practice, using numerical modeling. A consequence 
of this additional static potential is that cOz is different: 



-of 



MQV^zi 



(5) 



To obtain the constants ri, a, and zi, we use the Charged 
Particle Optics (CPO) numerical modeling software package 
to model the trapping potentials. The lattice trap used for 
our experiments has a hole diameter of h = 1.14 mm and a 
spacing between the centers of the holes of d = 1.64 mm. A 
square section of the rf electrode measuring 10 lattice sites on 



each side was used for this modeling; for larger sections than 
this, the effect of adding additional sites on the potentials 
near the center was negligible. From a simulation of the 
trap, we obtain the value of the geometric factors: ri = 
3.1 ±0.1 mm, a = —4. Oil. 3, and zi = 19 mm for a top plate 
15 mm above the rf electrode. Errors arise from the nonlinear 
least-squares fit used to obtain ri and a from the (discrete) 
simulated potential. In Fig. |3| we compare the numerical 
potential for the lattice trap to the analytical potential from 
the multipole expansion, indicating that near the minimum 
of a given potential well the multipole expansion gives an 
accurate approximation to the simulated pseudopotential. 



III. ATOMIC ION EXPERIMENT 

We have observed stable confinement of ^^Sr^ ions in a 
6x6 lattice trap, and verified the model of the trap discussed 
in Sec. II by measuring the secular frequencies of the ions for 
one particular lattice site. The rf electrode is cut from a 
stainless steel mesh from Small Parts, Inc., Part No. PMX- 
045- A. It is mounted 1 mm above a grounded gold electrode 
on a ceramic pin grid array (CPGA) chip carrier (Fig. [4|. 
An additional planar electrode (the top plate) is mounted 1 
cm above the rf electrode, to help shield the ions from stray 
charges. Electrical connections to both the rf and ground 
electrodes were made using a UHV- compatible solder from 
Accu-Glass (part no. 110796). The vacuum chamber was 
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FIG. 3: Fit to Eq[T]of the CPO computed pseudopotential in the 
z direction at the center of one well. The fit yields ri = 3.1 ± 0.1 
mm, a = —4.0 ± 1.3, and zi = 19 mm. 
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FIG. 4: (a) Schematic of the cross-section of the trap assembly. 
The trap electrode is held above the CPGA center pad on top of 
two 1 mm thick glass slides, (b) Photograph of the trap mounted 
in the CPGA. Connections for RF and GND are shown, as are 
the optional control electrodes for the x and y directions. 
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FIG. 5: (a) A cloud of ions (circled) intersects the detection lasers 
traversing the trapping region. The bright spots beneath the ions 
are laser scatter, (b) An ion crystal (circled) with a lifetime of 
0(15 minutes) is observed in the trap. 



baked out to a base pressure of 7xl0~^^ torr. 

The trap is loaded with ^^Sr^ by photoionizing a beam 
of neutral strontium produced by a resistive oven. This is 
a two-photon process at 460 nm and 405 nm that has been 
discussed in [29 . Resonance fluorescence is imaged onto a 
CCD camera by simultaneously addressing the main 422 nm 
transition and the 1092 nm repumper. Typical laser powers 
used are 10 /iW of 422 and 50 /iW of 1092. Ions were observed 
as both clouds and crystals (Fig. [5|. The cloud lifetime is 
quite short (0(10 s)), but a small crystal has been kept in 
the trap, illuminated with cooling light, for up to 15 minutes. 
This short lifetime is attributable to the vacuum pressure. 

A typical voltage of V = 300 V at Q/27T = 7.7 MHz was 
applied to the rf electrode using a power amplifier and helical 
resonator. Numerical modeling of the resulting pseudopoten- 
tial yields secular frequency values of a;^/27r = 300 kHz and 
uJz/'^Tr = 600 kHz. In order to test the model, we measure 
both secular frequencies as functions of the applied rf volt- 
age V. We also compute a trap depth of 0.3 eV, which is 
the energy required for an ion at the potential minimum to 
escape. 

Secular frequencies were measured for one site near the 
center of the lattice using the standard method of applying 
a low-amplitude (~0.02 V) oscillating voltage to the top plate 
at the motional frequency of the ions. When each vibrational 
mode of the ions is stimulated, their heating causes measur- 
able drops in the fluorescence intensity. This experiment was 
performed and compared to the model for several values of 
the drive voltage (Fig. |6|. Agreement is very good; mea- 
sured data points differ from the predicted values by at most 
5%, an error that results mainly from the approximation of 
the trap electrodes as perfect two-dimensional conductors for 
simulation. 

Although other sites near the center were also loaded, sec- 
ular frequency measurements are presented here for only one 



site of the lattice. These experiments answer our questions 
regarding the ability to construct, operate, and accurately 
model a two-dimensional lattice ion trap: our agreement of 
less than 5% could be made better by refining the simulation 
methods. 



800 
700 
600 
500 
400 
300 
200 
10C 



■ vertical frequency (model) 
vertical frequency (fit) 

■ radial frequency (model) «» 
radial frequency (fit) 




150 



200 



250 300 
RF Amplitude [V] 



350 



400 



FIG. 6: Secular frequencies as a function of rf voltage for one site 
of the lattice trap. Circles represent data points, dotted lines rep- 
resent linear fits to the data, and the solid lines are the predicted 
values from the model. The upper (red online) data are values of 
(jJz/'^TT, and the lower (blue online) are values of Ur/'^TT. 



IV. MACROION EXPERIMENT 

Another important test of the applicability of this lattice 
design to quantum simulations is the strength of interactions 
between ions in different wells. The charge to mass ratio of 
the strontium ions is unsuitable for this measurement in a 
lattice of this {d = 1.64 mm) spacing. Instead, a lattice 
trap was loaded with "macroions," aminopolystyrene micro- 
spheres with diameter 0.44 /im (Spherotech Part No. AP- 
05-10). The charge-to-mass ratio Q/m of macroions used in 
the experiment leads to observable repulsions between ions 
in neighboring wells, although it takes on a relatively wide 
range of values due to the fact that Q/m is not the same 
for every macroion. The use of macroions is also experimen- 
tally much less demanding than atomic ion trapping, because 
UHV pressures and laser cooling are not required. In fact, 
ions can be trapped in atmospheric pressure more easily than 
under vacuum, since air damping of ion motion increases the 
range of parameters suitable for stable trapping [30l |3T1 132] . 

Fig. [S] is a diagram of the experiment, which is an adapta- 
tion of the experiment in ^30^. The main components of our 
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Lattice Trap Schematic 
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FIG. 7: Experimental setup for the macroion experiment (Sec- 
tion IV). (a) 3-D schematic of lattice trap setup, (b) The lattice 
RF plate, as mounted for the microsphere experiment. The hole 
diameter is 1.14 mm and the hole spacing is 1.67 mm. The trap 
is supported by a printed circuit board. 
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FIG. 8: Trapping apparatus. The lattice trap is inside a plastic 
chamber which can be pumped down to ~ 1 torr. Macroions are 
loaded via the electrospray and the 4-rod trap, which extends 
through one side of the chamber over the surface of the lattice 
trap. 



apparatus are the electrospray system and the 4-rod loading 
trap. To load the ions in the lattice trap, we perform a mod- 
ifcation of the method in [33 , skipping the washing step. We 
prepared a buffer solution of 5 mL pure acetic acid, 26 mL 
IM NaOH, and 5% suspension microsphere solution. The 
buffer solution reduces spread in macroion charge. We son- 
icated the solution for 10 minutes to mix the microspheres 
evenly in the solution, added 30 mL of methanol, and again 
sonicated for 10 minutes. 

Compressed air, at a pressure of between 3 and 5 Psi, 
forces the buffer solution first through a 0.45 jam filter and 
then to the electrospray system. Here a copper wire at a 
voltage of 4 kV is inserted in the tubing and ionizes the so- 



FIG. 9: Image from above of ions in the lattice. The dark holes 
are the holes in the rf electrode; the grounded plane is 1.4 mm 
beneath them. Single macroions appear as white dots that are 
levitated above the plane of the rf electrode and are illuminated 
by 532 nm laser radiation at 5 mW. White dots on the surface of 
the rf electrode are due to stray light scatter. The left image was 
taken at V = 300 V and Q/27v = 1200 Hz and the right image was 
taken at y = 300 V and Q/27r = 1960 Hz. In the left figure, in 
the top well, two ions are shown repelling each other in the same 
well. 



lution as it passes. The ionized solution travels through a 
thin electrospray tip directed at a perforated, grounded elec- 
trospray plate and a 4-rod Paul trap just behind the electro- 
spray plate. The electrospray tips were made from capillary 
tubes, which are heated and stretched to produce narrow 
openings of 75-125 jam. 

As the solution enters the 4-rod trap, the methanol evapo- 
rates and the charged microspheres break into small clusters, 
the macroions. The 4-rod trap is driven at the drive param- 
eters of the lattice trap and extends through the wall of a 
plastic chamber over the lattice trap. Inside the chamber, 
the 4-rod trap extends 0.75 cm over the lattice trap and the 
bottom rod of the 4-rod trap rests 1 mm above the ground 
plate. 

The lattice trap is supported by standings inside the cham- 
ber, which can be closed on all sides to block air currents and 
can also be sealed and pumped down to ~1 torr. Glass slides, 
which are coated with InTi02 so that one side is conductive, 
act as the top plate. They allow a top view of the trap, and 
are supported approximately 15 mm above the RF plate. 
The ions are then confined approximately 0.25 mm above 
the plane of the RF plate. An image of the ions in the trap 
is given in Fig. [9j 

Typical initial loading parameters for macroions were 
n/27r = 1000 Hz and V = 250 V. We also applied a DC 
voltage of U = — 10 V to the top plate to improve the trap 
depth. Before studying ion-ion repulsion, we estimate the 
Q/m of the macroions by measuring their secular frequen- 
cies {cOz)- To do this, we apply a low- amplitude tickle to the 
top plate and observe the resonances directly on a video cam- 
era as ions rapidly oscillate back and forth. A measurement 
of ujz vs. 17 is shown in Fig. [lO| Using Eq. |5) we fit these 
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FIG. 10: a;^ vs. Q for an isolated macroion at a drive voltage 
of y = 255 V. The data for V and 2.5 V come from the same ion. 



data to obtain a charge-to-mass ratio of 1.9 x 10~^e/amu. 

We measured the Coulomb interaction of ions in neigh- 
boring lattice sites for six pairs of ions. In each pair, we 
measured the offset of each ion from the center of the well, 
as shown in Fig. [9j Note that while taking data on separa- 
tion of two ions, we ascertained that wells adjacent to those 
containing the ions were all empty. The effect of a third ion 
in an adjacent well is significant. 

A simple model of the interaction of two ions across wells 
is given as follows. An ion is confined by a force — mcj^xi, 
where Xi is the ion offset from the center of the well. Since 
generally xi <C (i, where d is the lattice spacing, the ion is 
approximately repelled by a force QiQ2/4:7reoS(P . Here s is 
a screening factor and (5i,2 are the charges of the first and 
second ion, respectively. The screening factor 5 < 3, where 
s =3 for an ion sitting at a height .25 mm above an infinite 
conducting plane. 

Using the expression for uOr derived from Eq. |3) 



(6) 



When Qi is not equal to (52, then the confining forces are 
characterized by different ujr and the two ions have different 
offsets from equilibrium. The ratio in offsets, if the masses 
are comparable (mi ~ ^2), should be {Q2/Q1Y' We ob- 
served exactly such an asymmetry between the offsets of the 
two ions, where typically Q2/Q1 is between 1 and 5. There 
may be additional small asymmetries due to edge effects and 
the presence of the 4-rod trap as well as differences in charge. 
Fig. pT] shows the displacement of a pair of ions as the rf drive 
frequency is varied, for one experimental run. The spread in 
charge to mass ratios and accordingly unknown values of 
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FIG. 11: Ion displacement from the well center for two ions in 
neighboring wells, as a function of drive frequency. Due to charge 
asymmetry, the maximum displacements of the two ions differ by 
a factor of ten. The drive voltage is 350 V. The model breaks 
down for large displacements (high trap frequencies); fits only 
include data below Q./2ti = 2500 Hz. 



Q and m for each ion (as in Ref. [30]) does not permit us 
to compare the observed repulsion to a theoretical model. 
Nevertheless, we have answered our third question: ion-ion 
interaction in a mm-scale lattice trap is observable by the 
mutual Coulomb repulsion of the ions, albeit for a different 
system than the atomic ions that would be used for quantum 
simulation. 



V. SCALING BEHAVIOR 

The lattice trap discussed in this paper provides a fairly 
straightforward method for realizing a two-dimensional array 
of trapped ions. The question remains: how useful could this 
system be for quantum simulation of two-dimensional spin 
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models [21], and in order to do this, what ion- ion spacing will 
be required? The scheme of Porras and Cirac is based on a 
pushing laser that exerts a state-dependent force on trapped 
ions that are coupled by their Coulomb interaction. In the 
limit in which the Coulomb interaction is small compared to 
the trapping potential (which is the case for lattice traps), 
the coupling rate between the ions is given by 

7 ^ (7) 

STreom^d^u;^' 

where F is the magnitude of the state-dependent force, e 
is the charge of an electron, m is the mass of each ion, d 
is again the ion-ion distance, and co is the trap secular fre- 
quency. For the lattice trap, uj of Eq. [7|is cj^- ^ is assumed 
to be due to a tightly-focused laser beam, and arises from 
a spatially-dependent AC Stark shift. For a 5 W beam of 
532 nm radiation that is focused from 50 /im to 3.5 jam over a 
distance of d = 50 /im, in traps operating at a; = 27r-250 kHz, 
we calculate a J coupling of 1 kHz, which should be observ- 
able if the dominant decoherence time is greater than 2it / J . 
Similar values can be obtained by using less powerful lasers 
closer to the atomic resonance; we use the 532 nm beam as 
an example only because of the readily-available solid-state 
lasers at this wavelength. The motional decoherence rate 
expected in microfabricated surface-electrode traps becomes 
negligible relative to the internal state decoherence time if 
the trap is cooled to 6 K; rates for the former have been 
measured at as low as 5 quanta/s [34]. Internal state deco- 
herence times depend on the specific ion being used and also 
on classical controls, but coherence times as long as 10 s have 
been reported [35l [36] . 

Unfortunately, the scaling properties of lattice traps do not 
favor such a low secular frequency at small ion- ion spacings. 
To maintain trap stability, the q parameter {q = ^^2^2 ) 
must be held constant near 0.3 as d varies. However, the 
trap depth D oc qV cannot be allowed to decrease too much, 
since traps of depth below about 100 meV have proven to 
be difficult to load. We note also that ri scales roughly as 
d. Therefore, the voltage must remain as high as it is for 
large traps, and the drive frequency and secular frequency 
CO must increase as 1/d, since co (x qCl. According to Eq. ([7|), 
the increased trap frequency erases the gain of placing the 
ions closer together. Indeed, it appears in this regime that J 
actually increases linearly with a result noted in Ref. [27 . 

Greatly increasing the trap size is not only impractical, 
but may render the width of the ground state wave func- 
tion of each trapped ion comparable to the laser wavelength, 
leaving the system outside of the Lamb-Dicke confinement 
regime. While some gains might be made from using the 
stronger field gradients of a standing wave configuration for 
the 'pushing" laser, it is clear that the scaling of co with 1/d 



is a discouraging feature of lattice traps. Of course, other 
options should be explored; a pressing question is how to 
modify the lattice trap design to allow for low motional fre- 
quencies for the directions along which the ions interact. One 
simple possibility might be to decrease the drive voltage V 
(and consequently the trap depth) once the trap is loaded 
with ions and they have been laser-cooled to a temperature 
much lower than the trap depth. 



VI. CONCLUSIONS 

In conclusion, we have proposed a design for a layered 
planar rf lattice ion trap which contains many of the fea- 
tures desirable for quantum simulation, including the abil- 
ity to control the structure of the 2D lattice of ions. We 
mentioned above that in the present design the structure of 
the ion lattice can be controlled by the intersection of the 
Doppler cooling, photoionization, and atomic beams; how- 
ever, future realizations of the trap might also include indi- 
vidually controllable dc electrodes beneath each lattice site 
that could be used to eject unwanted ions. Also, in a future 
trap, the rf electrode could be specially fabricated to include 
only desired lattice sites. 

Of course, the square lattice used in this paper is only 
one possible geometry. Other lattices, including hexagonal 
ones, could also be used for quantum simulation. This would 
enable the possibility of observing spin frustration. As a first 
experiment, we envision trapping a triangular array of three 
ions and generating a spin-frustrated ground state. Although 
we have chosen to focus on the Porras- Cirac type spin model 
simulations, many of the schemes proposed for performing 
2-D quantum simulations could be realized in such a trap 

[HlSniEIlEaElEl]. 

Our implementation in a mm-scale lattice trap is the first 
demonstration of stable confinement of ions in such a trap, 
and meaurement of the secular frequencies has confirmed our 
theoretical models of the trap. These were essential questions 
to address before we can move forward towards traps of this 
type on the scale of tens of microns, at which point they could 
become useful for quantum simulation (provided the secular 
frequencies can be kept low enough, as discussed above). It is 
also crucial to be able to measure interactions between ions, 
which we have done using charged microspheres. We hope 
that this work will stimulate further efforts towards measur- 
ing interactions between atomic ions in a lattice ion trap, 
paving the way for two-dimensional quantum simulations. 

We gratefully acknowledge funding from the MIT Under- 
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